Characterizing the hydrodynamic response to the initial conditions 
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Abstract 

In hydrodynamics, the momentum distribution of particles at the end of the evolution is com- 
pletely determined by initial conditions. We study quantitatively to what extent anisotropic flow 
v„ is determined by predictors such as the initial eccentricity e„ in a set of realistic simulations, 
and we also show the importance of nonlinear terms in order to correctly predict V4. This knowl- 
edge will be important for making a more direct link between experimental observables and 
hydrodynamic initial conditions, the latter being poorly constrained at present. 



1. Introduction 

One of the most important probes of ultrarelativistic heavy-ion collisions is the anisotropic 
flow. Event-by-event hydrodynamics fit] provides a natural way of studying anisotropic flow 
and its fluctuations. In such calculations, one supplies a set of initial conditions (IC), evolves 
them through ideal or viscous hydrodynamics, and finally computes particle emission at the end. 
One can then write the azimuthal distribution of outgoing particles in a hydrodynamic event as a 
Fourier series 

— —^l+2Y^v„co^[n{(p-^>n)], (1) 

^ n 

or equivalently |e""^} = v„e""''^", where {■ ■ •} is the average in one event. Since the largest 
source of uncertainty in these calculations is the initial conditions, it is useful to identify which 
properties of the initial state determine each observable. This knowledge can then allow to 
constrain the initial state directly from data. 

It is well known that the initial average profile of non-central collisions is almond-shaped. 
One can obtain the eccentricity and the direction of the initial profile computing £2^'^*- - 
in the CM frame, where £2 is the participant eccentricity, is the participant 
plane, whose direction is the minor axis of the ellipse, and {• ■ ■} denotes an average over the ini- 
tial density profile. The hydrodynamic expansion produces an elliptic distribution of particles, 
i.e. £2 0^ V2, and the maximum of the distribution of particles is aligned with the direction of 
the steepest energy gradient in the IC, i.e 02 = *I'2. Nevertheless, due to the finite number of 
nucleons colliding in each nucleus-nucleus collisions, there are fluctuations in the initial profile, 
thus the relations V2 oc £2 and *I'2 - ^2 may no longer be valid in event-by-event hydrodynamics. 
However, the study of these relations for several models of IC in event-by-event hydrodynamics, 
showed that they are reasonably satisfied L2«^ A-tSJ. 

Preprint submitted to Nuclear Physics A November 1, 2012 



Symmetry considerations have been used to argue that higher-order flows (v„ with n > 3) 
should also be created by an anisotropy, e„. It has been shown that this is also true for n = 3, but 
it is not valid for n - 4,5 Nevertheless, Teaney and Yan 16|] have introduced a cumulant 

expansion in the initial density profile, in which s„ is only the first term in an infinite series, and 
they have suggested that the hydrodynamic response may be improved by adding higher-order 
or nonlinear terms. 

In order to understand which properties of the initial state determine v„ and 'I',,, we propose a 
simple quantitative measure of the correlation between (v„, ^F,,) and any proposed predictor, such 
as (e„, <1)A We then use this to find better estimators for the anisotropic flows. Details can be 
found in ill- 



2. Characterizing the hydrodynamic response 



In previous works |13|,U,|5|], the correlation of the anisotropic flow with the initial geometry 
has been studied by visually inspecting two types of plots representing separately magnitude 
and direction: the - O,, distribution, and a scatter plot of v„ versus e„. Instead, we carry 
out a global analysis, studying the correlation between the entire vector both simultaneously and 
quantitatively. For a given event in a centrality class, we write 

v„e'"*" = kE„e""^" + 6, (2) 

where is a proportionality constant, and S is the difference between the calculated flow and the 
proposed estimator, or the error in the estimate. The generalized eccentricity e„, used throughout 
this work, is: e„e"'*" = - |r"e'"<*| / {r"), which was proposed in Ref. |3i|, and later showed to be 
the lowest term in a cumulant expansion of the initial energy density i^. The best estimator is 
defined as the one that minimizes the mean square error (|fi^|), where (• • •) is the average over 
events in a centrality class. Thus, the best linear fit is achieved when 

(e„v„cos[nCI'„ - <D„)]) 

isf) • 

We define the quality of the estimator by 

Quality = A;-^^. (4) 



The closer Quality to 1, the better the response and the smaller the rms error (|£^|). Note that a 
negative Quality means that *P„ and O,, are anticorrelated. 

We now present results for Au-Au collisions at the top RHIC energy using the hydrodynamic 
code NeXSPheRIO 111] . NeXSPheRIO evolves initial conditions generated by the event generator 
NeXus, solves the equations of relativistic ideal hydrodynamics, then emits particles at the end 
of hydrodynamical evolution using a Monte-Carlo generator. For this work, we generated 150 
events in each of the 10% centrality classes, and we added 115 events in an extra class for 
events with zero impact parameter - particles from the pseudorapidity interval I77I < 1 were 
used. Eccentricities are obtained by averaging over the initial transverse energy density profile 
at z=0. Our hydrodynamical calculations contain fluctuating initial flow, as well as longitudinal 
fluctuations, thus the final flow measured is not entirely determined by the initial transverse 
geometry. In this sense, these results represent something of a worst-case scenario. 
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Fig. [U displays the Quality of the hydrodynamic response for elliptic and triangular flow as 
a function of centrality. As expected, the elliptic flow is driven to the almond-shaped overlap 
area for non-central collisions, where is approximately coincident with ^2- Even for central 
collisions, where V2 comes from fluctuations, the best estimator given by Eq. (|2]) is able to capture 
the physics of the elliptic flow fluctuations. The Quality for n = 3 is not as close to 1 as for the 
elliptic case, but one can see that the triangularity £3 is still a good estimator of the triangular 
flow. We have also tested the triangularity definition proposed by Roland and Alver |8] (not 
shown, see Ref. [TJ), where is replaced by in Eq. The triangularity with weight is a 
better predictor than with for most centralities. Since the weight gives more weight to outer 
shells of the initial energy profile, this means that the source of anisotropic flow moves inward 
as the collision becomes more peripheral. 
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Figure 1: Quality of the hydrodynamic response for elliptic (circle), triangular (square), quadrangular (triangle) and 
pentagonal flow (star) for different centrality bins. The leftmost points correspond to events with zero impact parameter. 



The estimators for n - A,5 differ qualitatively from our results for V2 and v'3, where £4 and 
£5 only give reasonable predictions for central collisions. Moving to peripheral collisions the 
Quality decreases and even become negative. This shows that e„ alone cannot be used to map 
the hydrodynamics response to the initial geometry for n = 4, 5. 



3. Finding better estimators 

One can hope to improve the estimator adding more terms in Eq. (|2]i, e.g. 



v„e 



- ksne'""^" + k'sy"^'" + &, (5) 



where e'„ and O^, are other quantities determined from the initial density profile (for instance, 
the next higher cumulant). The procedure to obtain the better estimator is similar: minimize the 
mean-square error (|£^|) with respect to k and k', take the average over events in each centrality 
class, then insert the values of k and k' in Eq. ^ and compute 

ilksne""^" + k' Ey"'^'"\^) 
Quality = . (6) 

One can show that the quality is higher than with just one term, Eq. (01). 

We want to apply the improved estimator to the quadrangular flow (for other flows see Ref. 
101), but we need to know which are £', and O',. Taking an elliptic initial condition, i.e an elliptic 
profile with only £2 4^ 0, then evolving it through ideal hydrodynamics, and computing the 
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particle emission, one will obtain, using Eq. ([T]l, the anisotropic flows and the event-plane angles. 
Since there is only £2 one should expect only v'2 + 0. The results for this initial condition 
are: V2 ?i as a response to and ^2 - <t2, = when n odd (by symmetry), but there is 
a quadrangular flow, v\ + 0, and *1'4 is aligned with ^2- Since there is no £4, so where does V4 
come from? The answer is, V4 is induced by the almond shape. 

As quadrangular flow can also come from £2, we want to investigate, first, v\ solely as a 
response to the almond shape. Thus, the natural estimator is v^e'^"^* - kie^e'^'^^'f- + £ - the 
square in the first term of rhs is used to preserve the rotational symmetry. The result can be seen 
in the Fig. |2] For mid-central collisions, where £2 is large, the nonlinear term is important, yet 
this estimator is not as good as previous estimators of V2 and V3. Therefore, to obtain a better 
estimator for v'4 is necessary to combine both contributions, thus using Eqs. (|5]l and (|6]l. Fig. |2] 
show that the resulting estimator is good for all centralities. For V5, the best estimator has linear 
e nonlinear terms: 65 and s^et, |7]. 
We have defined a quantitative measure of the quality of estimators of v„ from initial condi- 
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Figure 2: Quality for 3 different estimators of v\. Tlie combined estimator (square) is the best estimator. 

tions in event-by-event hydrodynamics. V2 and V3 are well predicted by £2 and £3, but for V4 
and V5 it is necessary nonlinear terms. These results provide an improved understanding of the 
hydrodynamic response to the initial state in realist heavy-ion collisions. 
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